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Abstract 

In supersymmetric (SUSY) models with the gravitino being the lightest SUSY 
particle (LSP), the SUSY breaking scale (i.e., the gravitino mass) could be deter- 
mined by measuring the lifetime of the next-to- lightest SUSY particle (NLSP). How- 
ever, for an ultralight gravitino of mass of 0(1) eV, which is favored cosmologically, 
the determination of the SUSY breaking scale, or the gravitino mass, is difficult 
because the NLSP decay length is too short to be measured directly. Recently we 
proposed a new determination of the gravitino mass by measuring a branching frac- 
tion of two decay modes of sleptons. In this paper, we investigate the prospects for 
determining the gravitino mass at LHC. For demonstration we take some explicit 
gauge-mediation models and show that the gravitino mass can be determined with 
an accuracy of a few 10% for an integrated luminosity 10 — 100 fb _1 . 



1 Introduction 



The presence of a gravitino is the most fundamental prediction in supergravity pQ and its 
mass m.3/2 is an important parameter to determine the supersymmetry (SUSY)-breaking 
scale. The gravitino mass is predicted in a wide-range region, m 3 / 2 = 0.1 eV — 100 TeV, 
depending on mediation mechanisms of the SUSY breaking to the SUSY standard-model 
(SSM) sector. The lowest mass region, m 3 / 2 = 0.1 eV — 10 eV, is very interesting in 
particular, since there is no astrophysical and cosmological gravitino problem at all in 
this mass region [2J. If it is the case, the gravitino is the lightest SUSY particle (LSP), 
and the next-to-lightest SUSY particle (NLSP) decays into the gravitino inside defectors 
if it is produced at collider experiments. Therefore, we may have a signal for the NLSP 
decay into the gravitino at future collider experiments such as LHC Furthermore, we 
may determine the gravitino mass from the mass and lifetime of the NLSP. However, for 
a light gravitino of mass m 3 / 2 ^ 10 eV the NLSP lifetime is very short; for instance, the 
decay length is given by ctnlsp — 0.55 /zm(m 3 / 2 /l eV) 2 (mNLSp/200 GeV) -5 for a slepton 
NLSP, which is difficult to measure at collider experiments. 

In a recent article [3], we proposed a new method to determine the gravitino mass 
(i.e., the SUSY-breaking scale) for a very light gravitino, by comparing the branching 
fractions of two decay modes of sleptons, t\ — > t\ + r + t and t\ — > i + G3/2, instead of 
using the lifetime of the NLSP. (In this proposal we assume that the NLSP is the lighter 
stau f\. £1 denotes the lighter smuon or selectron, £ is \i or e, and G 3 / 2 is the gravitino.) 
However, we did not examine if this method works at LHC. The purpose of this paper 
is to show that the above method to determine the gravitino mass is indeed effective at 
LHC for a certain parameter region of the SUSY-particle spectrum. To demonstrate our 
point, we adopt simple gauge-mediated SUSY breaking (GMSB) models pE] in this paper. 
However, our mechanism for the determination of the gravitino mass is applicable to any 
gauge-mediation model as long as sleptons are the lightest next to the gravitino LSP and 
m 3 / 2 < 10 eV. 
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2 Measurement of the branching fraction for the slep- 
ton decays 



In this section we explain the basic idea of our method. We assume that the NLSP is 
the lighter stau fi, and the lighter smuon and selectron (collectively denoted by £±) are 
heavier than the stau but lighter than the lightest neutralino Xi-> m fi < m i 1 < m x<{- ^ n 
this case, l\ have two dominant decay modes. One is the decay into the gravitino, and 
its decay rate is related to the gravitino mass: 

where Mp = 2.44 x 10 18 GeV is the reduced Planck mass. The other decay mode is 
the three-body decay into the lighter stau, i\ — > f 1 ± r T £. As pointed out in the previ- 
ous work PJ, if one can observe both of these decay modes, the gravitino mass can be 
determined. This is because the gravitino mass is written as 

^ /2 = /"r^bodyA — i — (2) 

3/2 487TMJ, VlVbodJ r 3 _body 1 ; 

and r 3 _bod y is calculable once relevant SUSY particles' masses are known. Thus, we may 
derive the gravitino mass m 3 / 2 by measuring the branching fraction r 3 _ body /T2_body 

Here we should comment on the measurement of the three-body decay rate r 3 _body 
In principle, r 3 _ body depends on various SUSY parameters. However, in most of GMSB 
models with a light gravitino (m 3 /2 = 0(1) eV), the lighter stau, smuon and selectron 
are approximately right handed sleptons, f\ ~ f R , l\ ~ £ R , and the lightest neutralino 
is almost bino, Xx ~ B. Therefore, we assume it is the case in the following discussion. 
Under those approximations, the three-body decay rate is given by [5] 

1 «EM /-l-(^+^) 2 



3-body - m Iit --^-J o dxx(l-x + r T -r f ){l-x + r s ) 



X\J (1 — x) 2 + r 2 . + r| — 2(1 — x)r 2 — 2(1 — x)r~ — 2r 2 r 2 
(l-a;) 2 (r|-l + x) 2 : 



(3) 



where r B = m^a/mi ,r? = m^jmi R ^r T = m T /mi R , and the lepton masses and m e 
have been neglected. Here, we have defined r 3 _bod y = T(£^ — > fi~T~£) + T(£ R — > ff r + £). 
We have checked that the approximation Eq.(j3]) is quite good and we can reproduce the 
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Figure 1: The contour plot of r 3 _ body . The left-side figure is for nig = 100 GeV, and the 
right-side is for nig = 200 GeV. The contour parameter is (0.1, 1, 10, 100) eV from left 
to right. 



true r3_bod y with an accuracy of factor 30% via the approximation Eq.([3]) for simple 
GMSB models explained in the next section. Therefore, by measuring the three masses 
m^o, rng , and m fl , one can estimate the r 3 _ body with a good accuracy. Contour plots 
of r 3 _bod y calculated by Eq.Q is shown in FigJTJ Comparing them with Eq.(JTJ), one can 
see that the two-body decay rates r 2 _bod y are comparable to the three-body decay rates 
r3_body in a certain parameter space. 

Let us consider how we can measure the branching fraction r 3 _body/r2-body at LHC. 
Comparing the three-body decay 

£ R -> £ + r(soft) +n -> £ + r(soft) + r + G 3/2 , (4) 

with the two-body decay 

£ R - £ + G 3/2 , (5) 

we can see that the three-body decay accompanies a hard tau. Thus, as r 3 _bod y /r2-bod y 
is larger, the number of taus produced in SUSY-like events becomes larger. Therefore, 
one may think that the branching fraction can be estimated by counting the excess of the 
number of hard taus in SUSY-like events. However, this method is troublesome because 
of the difficulty of tau-identification and the enormous backgrounds. 

Here, we propose an alternative experimental method to measure the branching frac- 
tion r 3 _bod y /r2-body At LHC, £r are likely to be produced through Xi' s decays. Thus, 
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Figure 2: The Mf+g- distribution. Here we set m fl = 190 GeV, = 200 GeV,m^ = 
300 GeV and m 3 / 2 = 1 eV. The Mg+g- distribution from the 2-body decay have 
an endpoint at M 2 max = Jrr? § - mj R = 224 GeV, and 3-body at M 3 max = 

Jm| - m \ R \Jt - i m T + ^fi) 2 /^| fl = 63 GeV. In this case, r 3 _ body /T 2 _ body = 1-23. 



we consider the two decay chains, Xi ~^ I ~^ ^ ^G 3 / 2 an d Xi ~^ i 9 £r ~> ^ ± ^ =F 'rf 1 . 
The dilepton invariant mass from the former chain has distribution with a sharp edge at 
Me+e- = \/ m |o — tti~ . On the other hand, for the latter case the dilepton mass distribu- 
tion has a peak which has an endpoint at M#+i- = ^J m ^,o — m| ^Jl — (m T + mf 1 ) 2 /m'j . 
An example of the dilepton mass distribution is shown FigfSJ From the ratio of these two 
peaks' areas, we estimate r 3 _ body /T 2 _ body . 

We should note that the latter endpoint is crucial to determine the mass difference 
Am = mi R — {m T +mf l ) 1 which is one of the most important parameters for the calculation 
of r 3 _ body (see Fig. [1]). We also note that there are little background for the signal by 
virtue of the flavor subtraction technique, e + e~ + — e ± y u T , where each dilepton 

represents the final-state one and the sign between each dilepton event defines if the 
event is added or subtracted when booked in the histogram. 



Table 1: GMSB parameters of the models. 
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M (TeV) 
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" grav 


SPS7 
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Modell 


40 


80 


3 


13 
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Model2 


40 


80 


3 


10 


+ 


5 



3 Gauge-mediation models 

We consider a simple gauge-mediation model, where a SUSY breaking field S couples to 
N pairs of messenger chiral superfields, ip and ip, which transform as 5 and 5* under 
the SU(5)gut'- W = kipipS. The S field develops a vacuum expectation value k(S) = 
M + 9 2 F 1 where M is the messenger mass. With these conditions the low-energy spectrum 
of the SUSY particles including the gravitino mass are determined by 6 parameters, 
A = F/M, M, N, tan/3, sgn(/x) = ±1, and C grav [4]. The gaugino masses are generated 
from loop diagrams of the messengers and are given by, at the one-loop level, 

m a = ^A(a=l,2,3), (6) 

47T 

where A = F/M and «i = 5«em/(3 cos 2 6w)- Scalar masses, at the two loop level, are 
given by 



2^A 2 E^j C a (i), (7) 

where C a (i) are Casimir invariants for the particle (pi (Ci(i) = 3Y 2 /5). Here, we have 
omitted the higher order terms in an expansion in F/M 2 . The above gaugino and scalar 
masses are given at the messenger scale, and the physical masses should be obtained 
by solving the renormalization group equations. Finally, the gravitino mass is given by 
m 3/2 = (C grav /VS)(F/M P ). 

For numerical analyses in the next section, we choose a few model points: one is the 
Snowmass point SPS7 [6], and we also take two other model points to demonstrate the 
dependence on the model parameters. The GMSB parameters of those models are shown 
in Table. [Q In Figj3j those model points are shown in (m 3 / 2 , tan/?) -plane, together with 
contour plots of r 3 _body/r2-body The full mass spectrum of the models are shown in 
Table. EE 

1 Here, the mass spectrum is calculated by ISAJET 7.67 jS], and used in the event generation in the 
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Figure 3: The contour plot of r 3 _body/r2-bod y in the (777,3/2, tan /3) plane. We set A = 
40 TeV, M = 80 TeV, N = 3, sgn(/j) = +. In this figure, we calculate the mass spectrum 
and the decay rates by the program SOFTSUSY [7] as in Ref. [3]. The marks represent 
the model points we discuss in the text. 

4 Determination of the gravitino mass at LHC 

In this section we show that the gravitino mass can indeed be determined at LHC, by 
taking SPS7, model 1 and 2 as examples. In all analyses we use an event generator 
HERWIG 6.5 P]. 

Let us first consider the parton-level signatures of signal events. In FigJH the distribu- 
tion of the opposite charge dilepton invariant mass Mu is shown for the model lj^| Here 
the flavor-subtraction in the final state, e + e~ + fi + fi~ — e ± fi T , is adopted. As expected, 
one can clearly see the two peaks from two- and three-body decays of the slepton. 

following section. 

2 It seems that the matrix element of the three-body decay Ir — > t^t^I is not implemented in Herwig 
(i.e., the matrix element is taken to be constant). Note that the end point of Mu distribution from the 
three-body decay and the branching fraction is unchanged even if one takes it into account, and therefore 
only the shape of the first peak is affected. One may wonder that the reduction factor Eq. ([5]) due to the 
soft-lcpton cut Pt > 6 GeV is affected by this approximation, but we have checked that the inclusion of 
the correct matrix element does not change the result much. 
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Table 2: Mass spectrum 



Particle 


SPS7 


Modell 


Model2 


Particle 


SPS7 


Modell 


Model2 


g 


952.3 


952.1 


952.4 


G3/2 


0.77 eV 


0.77 eV 


3.85 eV 


u L 


902.1 


902.1 


902.2 


UR 


872.8 


872.7 


872.9 


d L 


905.8 


905.8 


905.8 


d R 


871.3 


871.3 


871.4 


b 2 


876.1 


875.9 


875.5 


61 


863.8 


865.0 


866.9 


*2 


895.8 


896.1 


896.6 


t x 


811.1 


810.9 


810.2 


V 


251.1 


251.1 


251.1 


V T 


250.6 


250.7 


250.9 


e L 


267.5 


267.5 


267.5 


en 


129.1 


129.1 


129.0 




268.8 


268.5 


268.0 


fx 


122.3 


124.0 


126.0 


xi 


160.1 


160.0 


159.8 


xt 


274.9 


275.0 


275.6 


x§ 


325.0 


325.6 


327.6 


xl 


388.8 


389.3 


390.9 


xt 


269.4 


269.1 


268.9 


xt 


391.2 


392.0 


394.0 


h° 


113.6 


113.4 


112.8 


H° 


389.7 


393.5 


400.1 


A 


389.5 


393.2 


399.5 


H± 


397.8 


401.4 


407.5 



However, events are subject, at LHC, to various cuts from experimental constraints and 
from reducing the backgrounds, and hence the real event distribution may not necessarily 
follow the ideal one. In fact, in many cases, the reduction factor 

^ # of dileptons with the cuts 
# of dileptons without any cuts 

tends to be non- negligible. If the reduction factor differs between for two-body and for 
three-body decay events, the Ts-body/^-body is affected by the cuts. 
In fact, we require that 

• the dilepton mass is formed only if one of the two leptons has Pp > 20 GeV, \rj\ < 2.5 
and the other has Pt > 6 GeVj3 < 2.5, 

where Pt and r\ denotes the transverse momentum and the pseudorapidity, respectively. 
This requirement affects the ratio r 3 _ body /r2_body Because of the small difference be- 
tween and fni R , the lepton from the three-body decay is typically soft. Therefore, 
the first peak from the three-body decay shrinks by the Pp cuts. The distribution of M« 
after the cuts is also shown in FigJH by the hatched histogram. As can be seen, only the 
events of three-body decay is substantially reduced. 



3 We have checked that our method works even if a harder cut Pt > 10 GeV is taken. 



s 




20 40 60 80 100 120 140 

M„ (GeV) 

Figure 4: The distribution of Mu- As expected, we can see two peaks of endpoints at 
M max _ 21.3 GeV and M 2 max = 94.5 GeV. The hatched histogram shows the distribution 
after the cuts on the lepton Pt- The peak from the three-body decay is clearly shown in 
the region of M u ^ 20 GeV, and in 60 GeV < M u < 90 GeV one can see the characteristic 
shape of Mu distribution from two-body decay. Note also that in Mu ^ 60 GeV there is 
another tiny peak which comes from the decay chain of the next to lightest neutralino, 
X2 — * ^ ±j ^r ~~ * ^ ± ^ T r 1 r, which has an endpoint at ^m^o — m| J^- — (jnfi + rn T ) 2 / rrv~^ = 
54.7 GeV. 

We have found that there is a correlation between the ratio of two peaks' endpoints 
M^^/M™^ = ^Jl — (m^ + m T ) 2 /m| and the reduction factor R 3 of the three-body 
decay as shown in Figl5l Here we have adopted various parameter regions in the simple 
GMSB model described in Sec. EJ We will use this correlation to reproduce the true 
branching fraction. 

Now let us discuss experimental signatures for each models in turn. For a detec- 
tor simulation, we use a package AcerDET-1.0 [TO]. We select events by the following 
requirements; 

• at least four jets with P T > 25 GeV, where r-jets are excluded. 
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Figure 5: The correlation between M™ ax /M™ ax and R% for three-body decay (left) and 
M™ ax and R2 for two-body decay (right). We have adopted simple GMSB models with 
parameters (N, A, M, tan/3) = ({3, 4, 5}, 30 TeV, 60 TeV, 15) and =({3, 4, 5}, 40 TeV, 
80 TeV, 18). In each model, we further vary the value of and m| fl as free parameters, 
imposing m^o > rn~ lR + 10 GeV and 771/ > 100 GeV. 

• missing transverse momentum Pr,miss > 100 GeV. 

• M eff > 500 GeV, where 

4 

M eS = ^Tj+^miss. (9) 

jets(^r) 



• two leptons with Pp > 20 GeV and |r/| < 2.5. 

We have checked that the first three cuts reduce the two decay chains Xi ~^ ^ ~^ 
£ ± £ T G3/2 and Xi ~ > ^ £r ~^ £ P^rfi almost equally and hence the ratio of the number of 
dileptons from the two- and three-body decay chains remains almost unchanged. The cut 
of two high Pt leptons slightly change the ratio, but not very much|f| We then form the 
dilepton invariant mass for opposite charge leptons which satisfy Pp > 20 GeV, \rj\ < 2.5 
and Pt > 6 GeV, \rj\ < 2.5^] As discussed above, this last requirement affect the ratio of 
number of events, and hence the reduction factor R should be taken into account. 

4 In fact, when a pair of sparticle cascades in a SUSY event end up with two decay chains where one 
is xi —> tit and the other is —> IrI — > fir£(soft)f , this event is likely to be cut by this 
requirement of two high Pt leptons. This reduces the number of three-body decay events, and decreases 
the resultant gravitino mass by about 20%. 

5 In our analysis, we do not take into account of miss-identification of soft leptons. 
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Figure 6: The distribution of Mu for the model 1. The hatched histogram is standard 
model background mainly from t-i production 



4.1 The model 1 

In FigJfJ the flavor subtracted (e + e~ + — e ± /i =F ) distribution of the dilepton mass 

after the cuts is shown for the model 1. Here, we have assumed an integrated luminosity 
£ = 10 fb" 1 . We also show the main background from the t-i events. We fit the data 
over 70 GeV < Mu < 120 GeV via a line smeared with a Gaussian, and obtain 

M max = g4 A ± q 5 Ge y_ ( 1Q ) 

To find M™ ax , we fit the data with a function 

h(x) = a(x - M)9(-x + M) +bx + c } (11) 
over 19 GeV < M u < 30 GeV, where x = M u . Then we get 

M 3 max = 21.1 ±1 GeV. (12) 
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The estimation of the error is done by 'eye'. Then, we find that the number of dileptons 



for the two-body decay is N 2 = 2876 ±201, and N 3 = 958 ±51 for the three-body decay|_|H 
From Figj5j we estimate the reduction factor of the three-body decay is R 3 = 0.20 ±0.05, 
and those for the two-body decay R2 = 0.90 ±0.02. Here, the errors are systematic. Then, 
we fincj® 

= Ml = (i.50 ± 0.15) 1 (13) 

Suppose that we know mg = 129.1±0.5 GeV which will be determined from additional 
observations like rrijt and rriju distributions |12j . We can then calculate the approximate 
Ts-body from Trig 1 M™ ax and M™ ax , by using Eq. (JHJ)- In the present case, we obtain 

IVbody = 0.2lig;o? eV (the true value is r 3 _ body = 0.22 eV). (14) 

Combining Eq. (fl3"|) and ()14p we derive the gravitino mass from Eq.([2]) as 



This value is in a good agreement with the expected one m 3 / 2 = 0.77 eV. 
4.2 The model 2 

In FigJTl the flavor subtracted (e + e~ + — e ± /i =F ) distribution of the dilepton mass 

after the cuts is shown for the model 2, for an integrated luminosity C = 100 fb^ 1 . We 
also show the main background from t-t events. In this case, an integrated luminosity 
C = 100 fb _1 would be necessary to determine the gravitino mass. 

We fit the data over 70 GeV < M« < 120 GeV via a line smeared with a Gaussian 
and find M 2 max = 94.1 ± 0.2 GeV. Also, we fit over 9 GeV < M u < 19 GeV via h(x) 
in Eq.flTTl) and find M 3 max = 13.2 ± 3 GeV. Then, we estimate N 2 = 41627 ± 1066 and 
N 3 = 1561 ± 784. We can see R 3 = 0.04 ± 0.02, R 2 = 0.90 ± 0.02 from FigJS Then, we 
obtain 

r^_Aift_ (aM±0 W*W*\- (16) 



r 2 -bod, N 2 R 3 ' V0.907 V0.04 



6 Here, those numbers are obtained by extrapolating the fitted line. It may be modified at lower Mu 
region depending on the detector performance (cf. [11]), but such an effect is small. 

7 The estimation of 7V3 could be affected by the contamination of the events from the decay of the 
second lightest neutralino Xo> but its effect is negligible (cf. Fig. |4j) . 

8 Here, the correlation between N2 and N3 is neglected for simplicity, taking those parameters to be 
independent. 
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Figure 7: The distribution of Mu for the model 2. The hatched histogram is standard 
model background. 



Supposing m~ lR = 129.0 ± 0.5 GeV, we estimate 

IVbody = 0.006±o;oo5eV (the true value is r 3 _ body 
which leads to 

R2 \ 5 ( R3 



0.006 eV), 



m 3/2 = (2.3±f; 



4.5\ 
4) 



0.90 



0.04 



eV. 



(17) 



(18) 



The expected value is m 3 / 2 = 3.85 eV. 
4.3 The SPS7 

In FigJHl the flavor subtracted (e + e~+yU + /i~— e ± yU =F ) distribution of the dilepton mass after 
the cuts is shown for the SPS7 model point, for an integrated luminosity C = 100 fb _1 . 
We also show the main background from t-t events. 

We fit the data over 70 GeV < M« < 120 GeV via a line smeared with a Gaussian 
plus a peak from Z°-boson decays and find M™ ax = 93.5 ±1.4 GeV. Also we fit the data 
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Figure 8: The distribution of Mm at SPS7. The hatched histogram is standard model 
background. The peak around 90 GeV is due to the Z° decays. 



over 24 GeV < M u < 35 GeV via h(x) in Eq.fllTJ and find M 3 max = 25.7 ± 2 GeV. Then 
we estimate N 2 = 8860 ± 1281 and N 3 = 17265 ± 498. We can see R 3 = 0.35 ± 0.10, R 2 = 
0.90 ± 0.02 from FigjH Thus, we obtain 

r 3 -bod y _ N 3 R 2 , . ( R 2 \ ( R 



(5.01±0.84) — =- — - (19) 
r 2 -bod y N 2 R 3 1 ; \0.90y vO.357 1 ; 

Assuming that we know mi = 129.1 ± 0.5 GeV, we estimate 

r 3 -bod y = 1.16lJ;g|eV (the true value is r 3 _ body = 1.11 eV), (20) 
which leads to 



"•3/^(o.4iiS:;D(^j-(^)"%v. pi) 



The expected value is m 3 / 2 = 0.77 eV. 
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5 Conclusion and Discussion 



We have investigated the prospects for determining the mass of an ultralight gravitino 
at LHC, by measuring the branching fraction of two decay modes of sleptons. We have 
performed detailed analyses by taking some specific GMSB models, and demonstrated 
that the proposed method can indeed work at LHC. Although we have taken simple 
GMSB models, our method works independently of details of GMSB models, as far as the 
two decay modes of sleptons are seen and the relevant mass parameters are all known. 

So far in this paper, we have assumed that the missing particle is the gravitino LSP 
and discussed how to determine its mass by measuring the two decay modes of sleptons. 
Conversely, one may argue that, if the two characteristic peaks of two- and three-body 
decays of sleptons are simultaneously seen in the dilepton invariant mass distribution, they 
themselves suggest that the missing particle is not a neutralino but the gravitino LSP, 
and therefore the underlying model is a GMSB model. One can then perform analyses 
as presented in this paper and determine the gravitino mass, or equivalently the SUSY 
breaking scale, which will be one of the most important physics target after the discovery 
of SUSY. 
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